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Background: Statistical measures of chaos have long been used in the study of chaotic dynamics 
in the framework of the interacting boson model. The use of large number of bosons renders 
additional studies of chaos possible, that can provide a direct comparison with similar classical 
studies of chaos. 

Purpose: We intend to provide complete quantum chaotic dynamics at zero angular momentum 
in the vicinity of the arc of regularity and link the results of the study of chaos using statistical 
measures with those of the study of chaos using classical measures. 

Method: Statistical measures of chaos are applied on the spectrum and the transition intensities 
of 0^ states in the framework of the interacting boson model. 

Results: The energy dependence of chaos is provided for the first time using statistical measures 
of chaos. The position of the arc of regularity was also found to be stable in the limit of large boson 
numbers. 

Conclusions: The results of the study of chaos using statistical measures are consistent with 
previous studies using classical measures of chaos, as well as with studies using statistical measures 
of chaos, but for small number of bosons and states with angular momentum greater than 2. 

PACS numbers: 21.60.Fw, 21.60.Ev, 21.10.Re, 23.20.Js 


I. INTRODUCTION 

The Interacting Boson Model (IBM) [T], apart from 
being successful in describing the low-lying levels and 
electromagnetic transition intensities of even-even heavy 
nuclei, has also been used in studying transitions [a |3] be¬ 
tween the different dynamical symmetries of the model, 
by changing its parameters. The model is known to pos¬ 
sess 3 dynamical symmetries, namely U(5), 0(6) and 
SU(3). When a system possesses a certain dynamical 
symmetry, it is completely integrable. Away from these 
dynamical symmetries, one would expect chaos. How¬ 
ever, this is not always the case. A situation of great in¬ 
terest is the notion of quasidynamical symmetry (QDS), 
i.e., the approximate persistence of a symmetry, in spite 
of strong symmetry breaking interactions ms]. Symme¬ 
try breaking can also be seen as the transition of a sys¬ 
tem from regular dynamics (exhibited by the presence of 
a dynamical symmetry) to chaos [HI |TU] . 

The interplay between regular and chaotic behavior in 
the context of the Interacting Boson Model (IBM), has 
been extensively studied by Alhassid and Whelan [nHS] 
and other authors HU, using both classical and quantum 
measures of chaos. In their study, they had found inte- 
grability at the three dynamical symmetry limits of the 
symmetry triangle [T7] of the IBM, namely, U(5), 0(6) 
and SU(3), as well as at the 0(6)-U(5) side of the tri¬ 
angle, due to the 0(5) symmetry known |18j to underlie 
the 0(6)-U(5) line. Away from these integrable regions, 
one expected chaotic behavior. However, the study of the 
interior of the symmetry triangle of the IBM, brought to 
the surface a region of nearly regular behavior [miiiiiii], 
connecting the U(5) and SU(3) vertices, known as the 


“Alhassid-Whelan arc of regularity” (AW arc). 

The increased regularity observed in the region of the 
Alhassid-Whelan arc, as well as the locus of the arc, have 
been studied using several different techniques. 

1) Information entropy of the wavefunctions is a mea¬ 
sure that quantifies the eigenstate localization of a par¬ 
ticular IBM Hamiltonian in different symmetry bases as¬ 
sociated with dynamical symmetries, linking it with the 
degree of regularity nans]. It has indeed been found 
[l9l [20] that increased localization in the symmetry bases 
occurs on the Alhassid-Whelan arc of regularity. 

2) The line corresponding to the degeneracy of the Oj 
and 2^ states within the symmetry triangle of the IBM 
has been found to closely follow the arc of regularity 
mini]. Furthermore, more than twelve nuclei exhibiting 
this behavior were placed on the arc, providing an exper¬ 
imental confirmation of its existence [2I1I22]- Later on, 
the approximate degeneracy of the Oj and 2^ states has 
been related to the degeneracy of the /3 and 7 bandheads, 
its locus has been determined through the intrinsic state 
formalism, and has been found to be located very close 
to the Alhassid-Whelan arc [131 ni] • 

3) The dynamics of 0+ states have also been consid¬ 
ered, using the nearest neighbor spacing distribution of 
O'*' states, in order to demonstrate the semiregular na¬ 
ture of the arc of regularity, in agreement with results 
obtained using classical measures based on Poincare sec¬ 
tions [13]. In particular, a bunching pattern of 0+ states 
has been found [H] on the arc, similar to the bunching 
pattern seen along the 0(6)-U(5) side of the triangle [25] . 
which is known for its regular dynamics [18] . 

4) The line corresponding to the degeneracy of the 2 ^ 
and 2+ states has also been found [IB] to closely follow 
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the arc of regularity. In the large boson number limit this 
degeneracy also guarantees the degeneracies predicted by 
the SU(3) symmetry in the first few bands lying lowest 
in energy. This can be considered [SB] as a sign of an 
underlying SU(3) quasidynamical symmetry, however its 
validity is limited to low lying states and to the large 
boson number limit. 

5) A line within the symmetry triangle of the IBM, also 
lying close to the arc of regularity, has been obtained [ 2 S] 
using a contraction of the SU(3) algebra to the algebra 
of the rigid rotator. This finding is related to the ground 
state band alone and has again been obtained in the limit 
of large boson numbers. 

6 ) Recently, families of high-lying regular rotational 
bands have been found [HHZ] in the IBM framework, 
occurring even in nuclei far away from the SU(3) dynam¬ 
ical symmetry and leading to increased overall regularity. 

The present paper offers a closer view on the quantum 
chaotic dynamics at the vicinity of the arc of regulatiry, 
without attempting to elucidate the nature of the sym¬ 
metry underlying the arc. Taking advantage of the new 
code IBAR [30l [3T] , which can handle up to Nb = 1000 
bosons, the present study is focused on states with zero 
angular momentum (J = 0), but can be easily extended 
to other angular momentum values. The main objectives 
are described here. 

1) To determine any energy dependence of statistical 
measures of quantum chaos in O'*" states and in B{E0) 
transition strengths. This study is for the first time feasi¬ 
ble, due to the good statistics allowed by the large num¬ 
ber of bosons used. 

2) To examine the stability of the location of the arc 
of regularity within the symmetry trianlge of IBM with 
changing boson number. The IBAR code allows for the 
first time to examine this stability for large boson num¬ 
bers. The question of stability is of importance in rela¬ 
tion to the empirical evidence [ 2 T 1 [ 22 ] found for the arc of 
regularity, since different nuclei are described by different 
boson numbers, in the region of deformed nuclei already 
used Elllig the boson number being close to 14. 

In Section II, the IBM Hamiltonian is described. The 
fluctuation measures, used to study the quantum dynam¬ 
ics in the symmetry triangle of the IBM are introduced 
in Section III, while in Section IV the numerical results 
are presented. An 0(6) line is considered in Section V, 
while in Section VI the discussion of the results is given. 


II. IBM HAMILTONIAN AND SYMMETRY 
TRIANGLE 

In the context of the IBM, low lying states in nu¬ 
clei can be described in terms of a monopole boson, s, 
with angular momentum 0 and a quadrupole boson, d, 
with angular momentum 2. The 36 bilinear combinations 
(s^s, djjs, form a U( 6 ) spectrum generating 

algebra. The three dynamical symmetries of the model, 
U(5), SU(3) and 0(6), which correspond to vibrational. 


rotational, and y-unstable nuclei, respectively, are placed 
at the vertices of the symmetry triangle, shown in Figure 
which is the parameter space of the model. 

In what follows we use the IBM Hamiltonian [EKTll 

US], 


H {v,x) = c 


md + 


V- 1 
Nb 


Qx-Q 


X ’ 


( 1 ) 


where hd = dN d is the d boson number operator, Qx = 
{s^d + d'^s) -I- is the quadrupole operator, and 

Nb is the number of valence bosons. The parameters 
( 77 , x) Etre the coordinates of the triangle and serve for 
symmetry breaking. 77 ranges from 0 to I, and x ranges 
from 0 to ~ —1.32. By varying 77 and y, the three 
dynamical symmetries of the model can be reached. U(5) 
corresponds to ( 77 , y) = ( 1 , 0 ), 0 ( 6 ) to ( 77 , y) = ( 0 , 0 ) 
and SU(3) to ( 77 , y) = (0, ^f^)- Numerical calculations 
of energy levels and B{E0) transition rates have been 
performed using the code IBAR [301131] which can handle 
up to Vs =1000 bosons. 


III. ELUCTUATION MEASURES 


In this section, the fluctuation measures or statistics, 
which are applied to the spectrum and the transition 
intensities of 0+ states, are introduced. However, be¬ 
fore proceeding to the quantum chaotic analysis of the 
spectrum or of the transitions intensities, the eigenval¬ 
ues should be unfolded. The reason is that the Gaussian 
Orthogonal Ensemble (GOE) requires that the average 
level spacing, S, of a spectrum in the limit V —>■ 00 
should be constant, however, the ordered sequence of 
levels {El, E 2 , ■. ■ T Ejx) produced, for example from the 
IBM Hamiltonian, forms a spectrum in which the low 
energy levels have consistently larger spacings than the 
high energy ones. In order to be consistent with the 
GOE requirements, one needs to unfold the spectra, that 
is to say, modify the spectrum, so that the average level 
spacing, S, is constant. The first step of unfolding is per¬ 
formed by constructing a staircase function of the data. 
A staircase function is the number of levels found be¬ 
low some specific energy. Then, a low order polynomial 
N{E) is fitted to the staircase function [32]. The un¬ 
folded energies, called normalized energies, are defined as 
Ei = N{Ei). With this mapping, the average level spac¬ 
ing of the spectrum of the normalized energies, the un¬ 
folded spectrum, becomes constant and specifically, equal 
to one, {S) = I. The fluctuation measures, are then ap¬ 
plied to the unfolded spectrum, which might have con¬ 
stant average level spacing, however, the spacings still 
show strong fluctuations. 

The statistical measures used for the determination of 
the fluctuation properties of the unfolded spectrum are 
the nearest neighbor level spacing distribution P{S) [ 55 ] 
and the A 3 statistics of Dyson and Mehta [511 ISi] ■ For 
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the quantum statistical analysis of the transition intensi¬ 
ties, the distribution P{y) [36l[37], a distribution, with 
ly degrees of freedom is applied to the B{E0) transition 
intensities. 

P{S) is defined as the probability that two adjacent 
energies differ by an amount of S. First, the normalised 
spacings. Si = Ci+i — et, are calculated and placed into 
bins, so a histogram of normalised spacings is produced. 
Then, the Brody distribution [33] is fitted to the his¬ 
togram 


= Aa(l-ha;)S'“exp(-a5'i+‘^), (2) 

where A is a scaling factor and a = r[(2-|-a;)/(l-|-CLi)]^+“. 
The value of w is found by fitting Eq. Q to the data by 
least squares. The Brody distribution takes the form of 
Poisson statistics when w = 0, which characterise a reg¬ 
ular system, and of the Wigner distribution when w = 1, 
which corresponds to a chaotic system. For intermediate 
cases, larger oj values imply more chaos. 

Spectral rigidity, is a measure of the deviation 

of the staircase function from a straight line, used to 
measure long range correlations. It was introduced by 
Dyson and Mehta [331 |3S] , who defined the function 

A 3 (a,L) = ^minA,B [N{E) - {AE + B)f dE, (3) 

where the constants A and B will give the best local fit to 
N(E) in the interval a < E < a + L and L is the energy 
length of the interval. For a random Poisson spectrum 
(regular case), A 3 takes the form 

Ar(T) = (4) 


For the GOE (chaotic) case there is an approximate ex¬ 
pression, for large L, 

AP'®(L) = ^{logL - 0.0687). (5) 


The exact expression, good for all L, is also known [35] . 
In fact, in order to calculate the A 3 statistics, in terms 
of the normalised energies ei, e 2 ,..., e„, the function of 
A 3 (a,L) is used, as given in Eq. (22) in [33], 


A 3 (a, L) 



1 \ 3n 

16 


1 


1 ( -‘2,1 + l)h ) , 


( 6 ) 


where = e^ — (a-|-L/2) is the measure of the normalised 
energies with respect to the center of the energy interval 
(a, a + L). A 3 (L) is calculated in energy intervals of 

length L, which span the whole normalised spectrum, 
once for intervals starting at a = 0 and once for intervals 
starting at a = L/2. Then, A^(L) is found as the average 


over all A^{a^L). The form of the function that should 
be fitted on the data is [33] 

Al{L) = A^^^[qL) + AP{{l-q)L). (7) 

The value of q is again found from the fitting of Eq. 
0 to the data points. For (7 = 0 , the regular case 
is reached, A^{L), while for q = 1, the chaotic limit 
emerges, A^‘^®(L). For intermediate values of q, the be¬ 
havior of the system is closer to chaos as q is closer to 
1 . 

The last distribution, P{y), where y is the relevant 
transition intensity, e.g. B{E0), is constructed in such a 
way that P{y)dy is the probability of finding an intensity 
in the interval dy around y. After proper normalization 
of the transition strengths y (see [111111131]), their loga¬ 
rithms are assigned to bins and a histogram of normalized 
transition strengths is produced. Then, the interpolating 
function [Ml [33 


P,y{y) = A 


2<y> 


''1'^ 2/2 ^exp{—vy/2 < y >) 


r(l) 


, ( 8 ) 


where A is a factor added for scaling reasons, is fitted to 
the histogram, through a least squares fitting, in order 
to find the best value of v. However, since y is the loga¬ 
rithm of the normalized transition strengths, one should 
change the variable y of the intepolating function of Eq. 
Q to 2 : = logig(j/) and use for the fitting the form of 
the interpolating function after the change of variable. 
When = 1, the interpolating function reduces to the 
Porter-Thomas distribution [Mj, which is the GOE case 
(chaotic). 


PGOE{y) = /c, \ ^ ■^exp{-y/2 < y >). (9) 

^/2^T <y> y/y 


For small values of v regularity is expected. There is no 
formal expression for the regular case. 


0 ( 6 ) 



FIG. 1: The IBM symmetry triangle and the points in the 
vicinity of the arc, where calculations were performed. 
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IV. NUMERICAL RESULTS 
A. Quantum chaotic dynamics of 0’*' states 

Measures of chaotic dynamics were calculated at four 
different points in the IBM symmetry triangle, as is 
shown in Fig. namely on the SU(3) vertex, a point 
on the arc of regularity having parameters (rj, x) = 
(0.632, —0.803) and at two points with the same rj ly¬ 
ing off the arc at (??, x) = (0.632,-0.7) and (? 7 ,x) = 
(0.632,-1.1). The value = 0.632 corresponds to C = 
0.7 in a different parametrization [29], and was chosen in 
order to be in accordance with the value used in Ref. |28| . 
The region of coexistence begins at 77 = 0.8 . The points 
on the arc are given by the expression xiv) = ^"^ 77 — ^ , 
which was found by calculating the values of x versus 77 , 
where a (a measure of classical chaos) is minimized | 20 | . 

In order to have good statistics, Nb = 175 was used 
for the fluctuation measures P{S) and A 3 (L), produc¬ 
ing 2640 0+ states. For the P{y) fluctuation measure, 
Nb = 50 was used, producing 54,756 possible transition 
strengths, B{E0)s, between the 234 O'*' states. The rea¬ 
son for selecting Nb = 50 instead of Nb = 175, was that 
in the latter case more than 6 , 000,000 possible transi¬ 
tion intensities are produced between the 2640 0+ states, 
which renders the calculation impossible to run in terms 
of time, while for Nb = 50, the run time and the statis¬ 
tics are more than satisfactory. The allowed number of 
B{E0)s differs from point to point in the symmetry tri¬ 
angle of the IBM and this number gets larger, as the 
dynamics of the point get more chaotic. Figure]^ shows 
the results for the three statistical measures P{S), ^^{L) 
and P{y) applied on the three points illustrated in Fig. 
accompanied by the results on the SU(3) vertex, a point 
with regular dynamics. 

All three different measures of fluctuations show con¬ 
sistent results. The SU(3) vertex is the more regular 
point, followed by the point on the arc of regularity, 
which appears to also have regular behavior. Then come 
the points labelled by x = —0.7 and x = —1.1 which 
are indeed chaotic, with the x = —0.7 point being less 
chaotic than the x = —1.1 point. The fitting of the inter¬ 
polating function of the last statistical measure, P{y), on 
the SU(3) vertex is rather poor, which reflects the fact 
that the expression of Eq. (|^ does not reduce to the 
regular case for any value of u. 


B. Chaos in 0^ states as a function of energy 

The study of chaos as a function of energy, using the 
quantum spectrum, is for the first time possible, due to 
the large number of bosons used, which allows for good 
statistics. The spectrum of the produced states is divided 
into equal parts and each part is studied using the three 
statistical measures. For Nb = 175, 2640 0'*' states are 
produced. The spectrum is divided into 8 parts of 330 


TABLE I: Spectra of O’*" states obtained from the Hamilto¬ 
nian of Eq. 0 for Nb = 175 and for the ( 77 , x) parameter 
sets (0, =^) [SU(3)], (0.632, -0.803) [arc], (0.632, -0.7) and 
(0.632, —1.1). In all cases the first 0^ state is set at zero and 
all other energies are normalized to the energy of the second 
0^ state, i.e. to the energy of the first excited O”*" state. For 
each of the 8 intervals of 330 states each, into which the spec¬ 
trum is divided, the energy of the highest state in the interval 
is shown. 


state 

SU(3) 

arc 

X = -0.7 

X = -1.1 

T 

0 

0 

0 

0 

2 

1 

1 

1 

1 

330 

29.19 

26.84 

25.39 

28.64 

660 

38.14 

37.37 

36.52 

38.38 

990 

43.30 

47.02 

46.82 

47.32 

1320 

46.49 

56.16 

56.57 

55.74 

1650 

49.58 

64.90 

65.94 

63.78 

1980 

52.74 

73.42 

75.00 

71.54 

2310 

55.87 

81.66 

83.79 

79.16 

2640 

59.00 

92.74 

94.28 

91.17 


TABLE II: Numerical values of the parameter w of the Brody 
distribution Rj(S) of Eq. for the 8 parts of each spectrum 
for Nb = 175. The states 1-330 are labelled as interval 1, the 
states 2311-2640 are labelled as interval 8. The value of uj 
obtained from fitting the whole spectrum, coming from Fig. 
2, is labelled by “total”. 


interval SU(3) arc x = —0.7 x = —1.1 
P(S) CJ CJ CJ UJ 


1 

0.30 

0.26 

0.51 

0.54 

2 

0.17 

0.50 

0.67 

0.82 

3 

0.14 

0.31 

0.67 

1.02 

4 

0.00 

0.16 

0.39 

0.85 

5 

0.06 

0.10 

0.18 

0.54 

6 

0.10 

0.22 

0.22 

0.35 

7 

0.14 

0.04 

0.23 

0.18 

8 

0.32 

0.03 

0.10 

0.16 

total 

0.17 

0.18 

0.34 

0.53 

states and each part 

is studied separately for its chaotic 


dynamics. The highest energy in each interval is shown in 
Table m In all cases the first O'*' state is set at zero energy, 
while all other energies are normalized to the energy of 
the second 0 ''' state, i.e. to the first excited O’*" state, 
thus rendering the parameter c appearing in Eq. 0 
irrelevant. 

The spectra shown in Table are the original ones, 
before any unfolding is applied to them. The division 
of the spectrum in parts of equal number of states, is 
the same as the division of the spectrum in parts having 
equal energy difference, due to the normalization of the 
spectrum which has led to neighboring energies differing 






5 


P(S) 



s 



SU(3) 


X= -0.7 


arc 
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L 
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FIG. 2: (color online) Results for the four points shown in Fig. for Nb = 175 and J = 0 states, using the fluctuation 
measures P{S) of E g. p] ) [55| and A 3 {L) of Eg. 0 mm , and for Nb = 50 and J = 0 states, for the fluctuation measure 
P{y) of Eg. 0 | 361 137|y The black dotted line shows the regular limit, the red dashed line describes the chaotic limit, while 
the blue solid line is the fit to the distribution. It is evident from the values of to, q and v, that the SU(3) vertex and the point 
on the arc are more regular compared to the other two points above and below the arc. 


on the average by 1. Figures[^and|^illustrate the results 
obtained for the nearest spacing distribution, P(S), and 
the spectral rigidity measure, A 3 (L), respectively. The 
first energy interval, i.e. the part with the states 1-330, 
appears at the top of the columns, while the last energy 
interval, the part with the states 2311-2640, appears at 
the bottom. The numerical results of w and q are dis¬ 
played in Tables [n| and |III[ 

The degree of chaos is not uniform in energy. The 
low energy part of the spectrum (the first energy inter¬ 
val) is always less chaotic than the closest higher parts 
of the spectrum (the second and third energy intervals), 
where the motion becomes apparently chaotic. However, 
at higher energies chaos decreases significantly, the spec¬ 
trum becoming almost regular at its highest part, even 
for the most chaotic point, y = —1.1. This behavior is 
common at all three points located at the vicinity of the 
arc. However, the SU(3) point, which is of course regular, 
seems not to follow the behavior of the others. Indeed, 
chaoticity falls in the middle part of the spectrum, and 
rises in the highest part, displaying exactly the opposite 
properties, although on the average the SU(3) point is 
less chaotic than the others, as indicated by the lowest 


values of w, q associated with it in Fig. 2. Another obser¬ 
vation is that the point on the arc is always less chaotic, 
in all energy intervals, compared to the points above and 
below the arc and that chaotic behavior is confined to 
fewer intervals at the arc or the SU(3) vertex. 

These results are in accordance with several previous 
studies of both classical chaos and quantum chaos; for the 
latter non-statistical measures of chaos have been used so 
far. A few examples are listed here. 

1) The energy dependence of regularity in the classical 
IBM has been studied for r] = 0.5 in Ref. [3S]. Increased 
regularity has been found at low energies and again at 
high energies, while reduced regularity has been observed 
in between. 

2) The energy dependence of regularity for a quantum 
IBM Hamiltonian has been studied, again for ry = 0.5, in 
Ref. [21]. Increased regularity has been observed in the 
AW arc region, both at low and at high energies. In this 
quantum study a visual method originally proposed by 
Peres [40] has been used, enabling a qualitative distinc¬ 
tion between regular and chaotic motion. 

3) The energy dependence of regularity in a quantum 
IBM Hamiltonian has been studied in a region corre- 
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P(S) P(S) P(S) P(S) 



FIG. 3: (color online) Results obtained for the nearest spacing distribution P(S) of Eq. ([^ [33] for each of the cases shown in 
Fig. 0 for Nb = 175 and, when the set of J = 0 states is divided into 8 sets of 330 states each. The states 1-330 (labelled as 
interval 1 in Table 0 appear on the top, the states 2311-2640 (labelled as interval 8 in Table 0 appear at the bottom. The 
black dotted line shows the regular limit, the red dashed line describes the chaotic limit, while the blue solid line is the fit to 
the distribution. 


spending to axially deformed ground states in Ref. |26|. Peres |40| has been used. 
Increased regularity, with strong occurence of SU(3)-like 
rotational bands, has been found at the AW arc at low 
energies and again at high energies, with reduced regu¬ 
larity occurring in between. Again the visual method of 
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L L L 

SU(3) x=-0-7 arc 


L 

z=-i.i 


FIG. 4: (color online) Results obtained for the spectral rigidity measure A'^{L) of Eq. 0 |34[ I35j , for each of the cases shown 
in Fig. [r for Nb = 175 and, when the set of J = 0 states is divided into 8 sets of 330 states each. The states 1-330 (labelled 
as interval 1 in Table III I appear on the top, the states 2311-2640 (labelled as interval 8 in Table III I appear at the bottom. 
The black dotted line shows the regular limit, the red dashed line describes the chaotic limit, while the blue solid line is the fit 
to the distribution. 


C. Chaos in 0’'’ states as a function of the 
parameter x 

The results of the quantum statistical parameters, w 
and q, as a function of the parameter x, determined at 14 


different points along the rj = 0.632 line of the triangle, 
are seen in Figurej^ First, being on the SU(3)-U(5) line 
of the triangle ((? 7 ,x) = (0.632, —1.32)) chaotic behavior 
is displayed, which keeps diminishing as one reaches the 
arc of regularity ((» 7 , x) = (0.632,-0.803)), where there 
































































TABLE III: Numerical values of the parameter q of the A'^{L) 
distribution of Eq. Q |34l I35| , for the 8 parts of each spec¬ 
trum for Nb = 175. The states 1-330 are labelled as interval 
1, the states 2311-2640 are labelled as interval 8. The value 
of q obtained from fitting the whole spectrum, coming from 
Fig. 2, is labelled by “total”. 


interval 

A3(L) 

SU(3) 

q 

arc 

Q 

X = -0.7 

q 

x = -i- 
q 

1 

0.01 

0.43 

0.67 

0.86 

2 

0.48 

0.84 

0.91 

0.92 

3 

0.52 

0.57 

0.79 

1.04 

4 

0.26 

0.32 

0.64 

0.97 

5 

0.40 

0.27 

0.37 

0.77 

6 

0.51 

0.41 

0.53 

0.60 

7 

0.54 

0.002 

0.53 

0.45 

8 

0.66 

0.26 

0.46 

0.44 

total 

0.39 

0.44 

0.63 

0.76 


X=-0.7 


arc 
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FIG. 6: Results obtained for the quantum statistical param¬ 
eters uj and q, for various values of at the three points in 
the vicinity of the arc of regularity. 
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-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

X 

FIG. 5: Results obtained for the quantum statistical param¬ 
eters UJ and q, for 14 different values of x along the rj = 0.632 
line of the triangle. The calculations were performed for 
Nb = 175 and J = 0. 


is a minimum. Then, as one moves to larger values of Xj 
chaoticity emerges again, but once more gives its place to 
regular behavior as one reaches the 0(6)-U(5) line of the 
triangle = (0.632,0)), where the 0(5) symmetry 

causes integrability [15] . 

The results for the O'*' states, are in complete agree¬ 
ment with the original work of Alhassid and Whelan, 
who found the existence of the arc of regularity using 
Nb = 25 and J > 2 angular momentum. A strong peak 
on the arc of regularity has also been noticed in Ref. [23] , 
where the authors studied the dependence on x at r]=0.5 
for the ratio of the number of regular trajectories to the 
total number of trajectories. 


D. Chaos in 0'*^ states as a function of the number 
of bosons Nb 

The results of the quantum statistical parameters, lu 
and g, as a function of the number of bosons Nb are seen 
in Figure In general, there is a drop of the values of 
the measured quantum statistical parameters w and q as 
Nb increases. For small number of bosons (till about 
Nb = 100) this drop is steep, while for larger Nb (larger 
than Nb = 175), the drop is very small and the quan¬ 
tum statistical parameters seem to have reached steady 
values. A reason for this steep drop, for small values of 
bosons, can be explained in terms of the total number 
of states for J = 0, for different number of bosons, seen 
in Table jTV] which affects the statistical results. For ex¬ 
ample, for Nb = 25 there are only 65 states, which give 
barely sufficient statistics, for Nb = 50 there are 234 
states which give better statistics, while for Nb = 175 
there are 2640 states which give more validity to the sta¬ 
tistical analysis of the eigenvalues. As the size of the 
system becomes larger, in the limit Nb ^ oo, uj and q 
seem to approach fixed asymptotic values. 

TABLE IV: Total number of states for J = 0, for different 
number of bosons Nb ■ 


Nb 

Number of states 

Nb 

Number of states 

25 

65 

175 

2640 

50 

234 

200 

3434 

75 

507 

225 

4332 

100 

884 

255 

5334 

125 

1365 

275 

6440 

150 

1951 

300 

7651 
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Despite the drop of the values of the parameters, as 
a function of Nb, the point on the arc has the smallest 
values of oj and q, for all Nb > 50, compared to the other 
two points, X = —0.7 and x = —l-lj the point x = —0.7 
being always less chaotic than x = —1.1 in the same Nb 
region. The situation is different for Nb = 25, where 
the value of w is greater at the position corresponding 
to the arc of regularity, than at the neighboring points 
(x = —0.7 and x = —1.1). This comes as a surprise, since 
in their work, Alhassid and Whelan had used Nb = 25 
and J = 2 or J = 10 in order to locate the arc of regular¬ 
ity. Probably, the failure to locate the arc, using the fluc¬ 
tuation measure P(S), in our case, has to be attributed to 
the marginal badness of the statistics, since for Nb = 25 
and J = 0 there are 65 states, while for Nb = 25 and 
J = 2 or J = 10, there are 117 and 211 states respec¬ 
tively. However, as the number of bosons increases, the 
limiting values are quickly reached and the position of 
the arc of regularity becomes stable. One should recall 
at this point, that the nuclei found to lie close to the 
arc of regularity [an [22] have boson numbers close to 
Nb = 14, therefore the location of the arc for these nu¬ 
clei might be different from the one corresponding to the 
Nb —>■ 00 limit. 

TABLE V: Numerical values of the parameter v of the dis¬ 
tribution Pv{y) of Eq. (|§ [3S1I3ZI, for the 8 parts of each 
spectrum for Nb = 50. The states 1-30 are labelled as inter¬ 
val 1, the states 211-234 are labelled as interval 8. The value 
of V obtained from fitting the whole spectrum, coming from 
Fig. 2, is labelled by “total”. 


interval 

P(y) 

arc 

V 

X = -0.7 

P 

x = -i- 

V 

1 

0.40 

0.62 

0.44 

2 

0.48 

0.81 

0.84 

3 

0.39 

0.60 

0.87 

4 

0.32 

0.49 

0.77 

5 

0.29 

0.35 

0.54 

6 

0.23 

0.26 

0.37 

7 

0.21 

0.21 

0.30 

8 

0.21 

0.27 

0.22 

total 

0.34 

0.42 

0.49 


E. Chaos in B{E0) intensities as a function of 
energy 

In the following subsections, a quantum chaotic study 
is carried out for the transition intensities between the 
0+ states. The chaotic behavior of B{E0) intensities 
as a function of energy is first presented. The 234 O'*' 
states occuring for Nb = 50 were divided into 7 parts 
of 30 states each and one part of 24 states. For each in¬ 
terval of 30 0+ states, the statistical measure P{y) was 


applied to the B{E0) intensities, produced by these 30 
O'*' states. The results are displayed in Table and in 
Figure The first energy interval, i.e. the part with 
the states 1-30, appears at the top of the columns, while 
the last energy interval, the part with the states 211 - 
234, appears at the bottom. The SU(3) point is missing 
from Table jVj for the following reason. As already men¬ 
tioned, for points possessing some dynamical symmetry 
and thus characterized by regularity, as the SU(3) point, 
the interpolating function of Eq. ^ is poorly fitted, 
since the number of B{E0) intensities occurring in this 
case is much smaller than the number of B{E0) inten¬ 
sities obtained at other points, characterized by chaotic 
dynamics. In addition, these relatively few B{E0) inten¬ 
sities are greatly dispersed. As a consequence, the study 
of chaos as a function of energy in the SU(3) limit, using 
Nb = 50, was statistically impossible. 

Again, the degree of chaos is not uniform in energy. 
The low energy part of the spectrum (the first energy in¬ 
terval) is always less chaotic than the next energy inter¬ 
val, where the motion becomes apparently chaotic. How¬ 
ever, at higher energies chaos decreases significantly, the 
spectrum becoming almost regular at the highest part of 
the spectrum, even for the most chaotic point, x = — 1 . 1 . 
This behavior is common at all three points located in the 
vicinity of the arc. Again, the point on the arc displays 
less chaoticity, in all energy intervals, compared to the 
points above and below the arc and chaotic behavior is 
confined to fewer intervals at the arc. What is fascinating 
though, is that the study of quantum chaotic dynamics, 
both on the transition intensities and the spectrum of 0 *^ 
states, show consistency in the degree of chaos as a func¬ 
tion of energy. In all cases, there is a jump in the degree 
of chaoticity as one passes the low part of the spectum, 
which is replaced by regularity as one gets at the upper 
part of the spectrum. 


F. Chaos in B{E0) intensities as a fnnction of the 
parameter x 

The results of the quantum statistical parameter v, as 
a function of the parameter x, along the 77 = 0.632 line of 
the triangle for Nb = 175 are seen in Figure]^ For these 
calculations Nb = 50 was used, which produced 54,756 
possible B{E0)s between the 234 0+ states. 

The results show the same general behavior, as those 
of the quantum statistical parameters uj and q. There 
is a characteristic minimum of the curve, at x = — 0 . 8 , 
revealing the minimum chaoticity of the arc. Chaotic be¬ 
havior is encountered at the regions below and above the 
arc of regularity, while at the 0(6)-U(5) line of the tri¬ 
angle (( 77 , x) = (0.632,0)), regularity emerges again. The 
study of chaoticity of the B{E0) intensities distribution 
gives results which are in complete agreement with the 
study of chaoticity of the B{E2) intensities distribution 
of the original work of Alhassid and Whelan. 
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P(y) P(y) P(y) 





y y y 

X=-0.7 arc X=-l-l 

FIG. 7: (color online) Results obtained for the statistical measure P{y) of Eq. ([^ |36[ I37j . for each of the cases shown in Fig. 

for Nb = 50 and, when the set of 234 J = 0 states is divided into 7 sets of 30 states and one set of 24 states. The states 
1-30 (labelled as interval 1 in Table |v| appear on the top, the states 211-234 (labelled as interval 8 in Table [v| appear at the 
bottom. The blue solid line is the fit to the distribution. 


V. STUDY OF CHAOS ON A LINE BASED ON 
0(6) PDS 

In addition to the notion of quasidynamical symme¬ 
try (QDS), i.e., the approximate persistence of a symme¬ 
try, in spite of strong symmetry breaking interactions Il¬ 


ls], the notion of partial dynamical symmetry (PDS) has 
been introduced including three different cases. 

In type I PDS part of the states possess all the dynami¬ 
cal symmetry, in type II PDS all the states possess part 
of the dynamical symmetry, while in type III part of the 
states possess part of the dynamical symmetry jS] . The 
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1.0 
0.8 
V 0.6 
0.4 
0.2 

-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

X 

FIG. 8: Results obtained for the quantum statistical param¬ 
eter V, for 14 different values of x along the r] = 0.632 line of 
the triangle. The calculations were performed for Nb = 50 
and J = 0. 


TABLE VI: Numerical values for Nb = 175 of the parameter 
uj of the distribution Pui{S) of Eq. (§ m, for the 7 points 
around the intersection (labelled as “on arc, on 0(6)”) of the 
0(6) PDS line with the arc of regularity. 



X 

V 

UJ 

below arc, on 0(6) 

-1.10 

0.60 

0.49 

on arc, on 0(6) 

-0.88 

0.54 

0.16 

on arc, right of 0(6) 

-0.98 

0.41 

0.23 

on arc, left of 0(6) 

-0.80 

0.63 

0.18 

above arc, on 0(6) 

-0.50 

0.38 

0.63 

above arc, right of 0(6) 

-0.50 

0.17 

0.55 

above arc, left of 0(6) 

-0.50 

0.57 

0.80 

linkage between QDS and 

PDS has 

been only 

recently 


clarified, by proving that coherent mixing of one symme¬ 
try (QDS) can lead to partial conservation of a different, 
incompatible symmetry (PDS) [45j . 

In Ref. [15], by using a measure of a fluctuations 
in a state dl, called Act$, a valley of almost vanishing 
Aag,s. fluctuations has been found, i.e., a valley where 
the ground state wave functions have a high degree of pu¬ 
rity with respect to the a quantum number of 0(6), pro¬ 
viding an example of an 0(6) approximate PDS of type 
III. This line begins from the 0(6) vertex and reaches the 
U(5)-SU(3) line of the triangle, intersecting, at a certain 
point, with the arc of regularity. 

As a preliminary test, we applied the nearest neighbor 
spacing distribution, on 0+ states, for Nb = 175 bosons, 
at various points on and around the 0(6) PDS line, in 
order to see whether the 0(6) PDS line induces regular 
dynamics to the O’*" states considered here. The results 
for the relevant parameter, ui, are presented in Table [VT| 
The only regular points are those located on the arc of 
regularity, including the intersection of the 0(6) PDS line 
with the arc. Therefore, the 0(6) PDS does not seem to 
induce regularity to the 0+ states considered here, a fact 


not unexpected, since the 0(6) PDS regards the ground 
state band, while the present study focuses on all 0+ 
states, the two sets having only one state in common. 

It is expected that in general PDS can lead to sup¬ 
pression of chaos, however this is expected to depend 
on the number of states exhibiting the PDS. In [35] a 
model Hamiltonian with an SU(3) PDS has been used, 
in order to examine whether there is suppression of chaos 
on the point of the PDS. Classical and quantum mea¬ 
sures of chaos have been used for three different values 
of the angular momentum (J = 2,10, 25) and Nb = 25. 
Minimum chaoticity has been found around the point 
where the SU(3) PDS was supposed to exist, but the 
minimum diverged slightly from this point for the two 
lowest values of angular momentum, while it was more 
pronounced and closer to this point for the larger value 
of angular momentum (J = 25). The same behavior has 
been seen for the average SU(3) entropy of the eigen¬ 
states of the model Hamiltonian, which also showed a 
more pronounced and closer to the point of the SU(3) 
PDS minimum for J = 25. This behavior has been ex¬ 
plained in terms of the increase of the number of soluble 
states with increasing J. In contrast, in |47j . where a 
large number of states had an SU(2) PDS, the classical 
and quantum measures of chaos showed a minimum ex¬ 
actly at the point where the SU(2) PDS was supposed to 
exist. 


VI. CONCLUSIONS 

In this paper, we presented the application of statis¬ 
tical measures of chaos on the energy spectrum and the 
transition strengths of excited 0+ states, in the context 
of the interacting boson model and more precisely on a 
line connecting a point on the U(5)-SU(3) leg of the sym¬ 
metry triangle with a point on the U(5)-0(6) leg. While 
the statistical measures of chaos reveal regularity on the 
U(5)-0(6) leg, due to the 0(5) symmetry known [18] to 
underlie this leg, the point on the U(5)-SU(3) leg seems 
to be the most chaotic point on this connection line. Be¬ 
sides the U(5)-0(6) leg, one more minimum is revealed 
on the transition line, located on the arc of regularity, 
which is in agreement with previous studies performed 
with smaller number of bosons and levels with angular 
momentum J > 2. 

Concerning the objectives posed at the end of the in¬ 
troduction, the following comments are in place. 

1) The study of chaos as a function of energy gave 
some intriguing results. This study has already been 
performed using both classical |23j and quantum mea¬ 
sures of chaos [24l [26], but never using the fluctuation 
measures of the quantum spectrum or the transition in¬ 
tensities, due to the bad statistics imposed by the small 
number of bosons used. The degree of chaos differs as 
the energy changes. The results, both for the spectrum, 
as well as the transition intensities, show that the low 
energy part of the spectrum has always more regular dy- 
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namics compared to the immediately higher energy parts, 
which are almost completely chaotic. However, as the 
energy increases, chaoticity decreases and regularity pre¬ 
vails at the highest parts of the spectrum. This is the 
general behavior for points studied in the vicinity of the 
arc. However, the point on the arc is always more regu¬ 
lar, compared to the other points, in all energy intervals, 
giving a confirmation for its regular character. 

2) Quantum chaos also depends on the number of 
bosons used. As the number of bosons increases, the 
degree of chaos converges to a steady value. Beyond 
Nb = 50, the relative chaoticity of different points is 
the same, for all numbers of bosons used in the context 
of this work. For example, chaoticity on the arc points 
is less than on the neighboring points, regardless of the 
number of bosons beyond Nb = 50, giving for the arc 
a steady position in the symmetry triangle as a function 
of boson number beyond Nb = 50. Below Nb = 50, 
however, the location of the arc appears to be sensitive 
in As, this fact affecting the efforts of finding nuclei ly¬ 
ing on or near the arc [HIES], since the boson numbers 
corresponding to these nuclei are around Nb = 14. 

While the fluctuation measures of chaos reveal semireg¬ 
ularity on the arc, they cannot reveal the nature of the 
approximate symmetry underlying the Alhassid-Whelan 
arc of regularity, a question which has been addressed in 
several recent studies [Ml ESI ESI ES] ) but still remains 
open. 

Recently, a line based on an 0(6) PDS was found in the 
symmetry triangle, extending from the 0(6) vertex to a 
point on the U(5)-SU(3) leg of the triangle, intersecting 


with the arc of regularity. Calculations of chaos around 
the intersection point show that the 0(6) PDS does not 
contribute to the development of regular dynamics for 
the 0+ states, a reasonable result since the 0(6) PDS 
under consideration regards only the ground state band 
of the nuclear spectrum. 

The extension of the present study to states with non¬ 
zero angular momentum is desirable. In this case, at¬ 
tention should be paid to handling degeneracies in the 
spectrum |15j . 

Study of the interplay of order and chaos has also been 
extended to the region of shape-phase transitions in the 
triangle. Shape phase transitions between different dy¬ 
namical symmetries, as well as critical point symmeties 
appearing at the relevant transition points, have been an 
active field of investigation over the last decade [HE]. In 
the framework of the interacting boson model, in particu¬ 
lar, a first order shape phase transition between spherical 
and deformed shapes is known to exist [48] , characterized 
by a phase coexistence region [33] ■ The appearance of de¬ 
generacies within the phase coexistence region [sniEi], 
as well as the evolution of order and chaos across the first 
order transition |52II55j . are emerging fields of investiga¬ 
tion. 
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